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We study the phase transition between the high temperature algebraic liquid phase and the low 
temperature ordered phase in several different types of locally constrained O(N) spin systems, using 
a unified constrained Ginzburg-Landau formalism. The models we will study include: 1, O(N) spin- 
ice model with cubic symmetry; 2, O(N) spin-ice model with easy-plane and easy-axis anisotropy; 
3, a novel O(N) "spin-plaquette" model, with a very different local constraint from the spin-ice. We 
calculate the renormalization group equations and critical exponents using a systematic e = 4 — d 
expansion with constant N, stable fixed points are found for large enough N. In the end we will 
also study the situation with softened constraints, the defects of the constraints will destroy the 
algebraic phase and play an important role at all the transitions. 

PACS numbers: 



I, INTRODUCTION 

It is well-known that sometimes constraints imposed 
locally can lead to stable phases with unusual long dis- 
tance correlations. For instance, at zero temperature the 
local constraints can give rise to quantum Bose liquid 
phase with gapless excitations analogous to photon or 
even graviton [1-5] . At finite temperature, in three di- 
mensional classical dimer model (CDM) and spin-ice sys- 
tems, the local geometric constraint grants the high tem- 
perature spin disordered phase an algebraic power-law 
correlation between dimer and spin densities (3-11 , in- 
stead of the standard short-range correlations according 
to the Ginzburg-Landau theory. This algebraic phase is 
usually called the Coulomb phase. In both the dimer 
models and spin-ices, physical quantities are defined on 
links of a bipartite lattice, and let us take the cubic lat- 
tice for simplicity. The ensemble of 3d CDM is all the 
configurations of dimer coverings, which are subject to 
a local constraint on every site: each site is connected 
to precisely m dimers with < m < 6 (denoted as 
CDM-m), and most studies have been focused on the 
case with m = 1. The ensemble of the spin-ice model is 
all the configurations of spins on the links, and the sum 
of the six spin vectors around each site is zero, which 
is usually called the ice- rule 12J. Although the spin is 
generically an 0(3) vector, due to the spin-orbit coupling, 
the spins in the spin-ice materials prefer to align along 
the links where the spins reside, therefore people usually 
treat the spins in spin-ice materials Ising spins [7] , hence 
the spin-ice model on the cubic lattice is mathematically 
equivalent to CDM-3. The power-law correlation of the 
Coulomb phase has been confirmed by numerical simu- 
lations and also neutron scattering in spin-ice ma- 
terials such as Ho 2 Ti 2 07 and Dy 2 Ti 2 0y [13 1. 

The partition function of the CDM-m is a sum of all 
the dimer configurations allowed by the constraint, with 
a Boltzmann weight that favors certain dimer configu- 
rations. To describe CDM-m concisely, we can intro- 



duce the "magnetic field" Bi^ — (n^ — m/6)rji with 
/i = +x, +y, +z, where r\i = (— l) 1 is a staggered sign 
distribution on the cubic lattice. The number m ^ is 
defined on each link (i, fx) between sites i and i + fi, 



and rii 



1, represents the presence and absence of 



dimer. Now the local constraint of the dimer system can 
be rewritten as a Gauss law constraint V • B = 0. The 
standard way to solve this Gauss law constraint is to in- 
troduce vector potential A defined on the unit plaquettes 
of the cubic lattice, and B = V x A 3, HI- Since vec- 



tor A is no longer subject to any constraint, it is usually 
assumed that at low energy the system can be described 
by a local field theory of A, for instance the low energy 
field theory of the Coulomb phase reads 
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which is invariant under gauge transformation A — > 
A + V/, / is an arbitrary function of space. In the 
Coulomb phase, the correlation of magnetic field B^ in 
the momentum space reads 



{B^{q)B v {-q)) ~ 6^ - 



(2) 



The spin-ice model can be straightforwardly general- 
ized to the O(N) case. We can define an O(N) spin 
vector S a with unit length J2 a (S a ) 2 = 1 on each link 
(i, jti) of the cubic lattice (Fig. Wp), with a = 1---N, 
and we assume that the largest term of the Hamiltonian 
imposes an ice-rule constraint [T3| for O(N) spins on the 
six links shared by every site: 
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The magnetic field formalism developed for the CDM can 
be naturally applied here, with 0? = %S? + , and the 
constraint Eq. [3] can be written as 



(4) 
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Then in the Coulomb phase the momentum space corre- 



lation between § a reads 



), 



(5) 



which after Fourier transformation leads to the same 1 /r 3 
power-law correlation as the dimer model. 

Another constrained spin system, which is very dif- 
ferent from the spin-ice is the classical plaquette model 
(CPM). In Ref. [lg, [l7[ a plaquette model engineered 



from the SU(4) Heisenberg model on the cubic lattice was 
studied, and one of the twelve unit square faces shared 
by each site of the cubic lattice is occupied by a pla- 
quette, which is physically the SU(4) singlet formed by 
four fermions at the corners. The exact SU(4) symmetry 
can be realized in cold atom system without fine-tuning 
Til 19| . The ensemble of the CPM-m is all the plaque- 
tte configurations on the cubic lattice, with exactly m 
plaquettes shared by each site (Fig. (2f)). To generalize 
this system to O(N), let us define O(N) spins S a on the 
unit faces instead of the links of the cubic lattice, and 
impose the constraint that the sum of spin vectors on all 
the twelve faces shared by every site be zero (Fig. 
When N = 1, this model is equivalent to CPM-6. In this 
case it is most convenient to introduce symmetric rank-2 
ViSiuv w i tn M> v = +&> +V> (h V* v ) de- 



tensor ' 



notes the unit square face shared by sites i, i+fJ., i + v and 
i + fi + v. In terms of <^„, the spin plaquette constraint 
can be written as 



0. 



(6) 



This spin-plaquette system also has an algebraic phase at 
finite temperature, with the following momentum space 
correlation: 



%<lv<lp<l<j 
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Besides the high temperature algebraic phases in the 
models discussed above, at low temperature spins are ex- 
pected to order according to the details of the Hamilto- 
nian. In this work we will study the phase transition be- 
tween the high temperature algebraic phase and the low 
temperature spin ordered phase in 0{N) spin-ice (section 
II), with both cubic symmetry and anisotropies, as well 
as O(N) spin-plaquette system (section III). In section 
IV we will study the situation with softened constraints. 
The nature of the transitions clearly depend on the low 
temperature spin order pattern, and in this work we will 
focus on one particular type of spin order, which has 
nonzero net (0°) and (</>„„) at large length scale. Our 
calculation will be based on e = 4 — d expansion. Due to 
the complexity of the calculation, in our current work we 
will keep the precision to the first order e expansion. 



II, O(N) SPIN-ICE MODEL 
A, Brief Review of the CDM 

Let us first give a brief review of the previously stud- 
ied phase transition of the CDM, or the spin-ice model 
with N = 1. The field theory Eq. [1] misses one im- 
portant piece of information: the magnetic field B and 
the vector potential A are both discrete. Therefore 
mathematically we should introduce "vertex operator" 
L v ~ J^fj, cos[2tt(A^ + a n)] to the field theory Eq. [TJ and 
a M is a nonzero background distribution of the vector 
potential, which is introduced for any nonzero m. The 
Coulomb phase is a phase where this vertex operator is 
irrelevant perturbativcly. The vertex operator will be- 
come nonperturbative and drive a phase transition when 
it is large, and in order to describe this phase transition 
one can introduce matter fields in the vertex operator 
which couple to the gauge field minimally: 



cos[2tt(V m 6> - A M -a M )], 



(8) 



9 is the phase angle of the matter field ip ~ e 10 . Due 
to the existence of the nonzero a^, the matter fields ip 
move on a nonzero background magnetic field, the band 
structure of the matter fields have multiple minima in 
the Brillouin zone, and the transformation between these 
minima encodes the information of the lattice symmetry 
[20| . Therefore in addition to manifesting the discrete na- 
ture of the gauge potential A^, the condensation of the 
matter field leads to lattice symmetry breaking, which 
corresponds to the crystal phase of the CDM. For in- 
stance the transition between the Coulomb and columnar 
crystal phases of the CDM-1 model is described by the 
CP(lj model with an enlarged SU(2) global symmetry 
0, [21j . This field theory is highly unconventional, in the 
sense that it is not formulated in terms of physical order 
parameters. It is expected that more general CDM-m 
models can also be described by similar Higgs transition, 
although the detailed lattice symmetry transformation 
for matter fields would depend on m. 

One might be tempted to describe the transition be- 
tween the Coulomb and columnar phases of CDM-1 
trough a Ginzburg-Landau approach. One can introduce 
an 0(3) vector with cubic symmetry anisotropy in fa- 
vor of six axial directions, and (y^) ~ ±1 represents six 
fold degenerate columnar order. However, the hedgehog 
monopole configuration of the 0(3) vector always in- 
volves a broken dimer i.e. a defect of the constraint. Then 
as long as we forbid the presence of the defects, this 0(3) 
model is monopole-free, and it is well-known that the 
monopole- free 0(3) nonlinear sigma model is equivalent 
to the CP(1) model (2^1, which has very different crit- 
ical exponents from the 0(3) Wilson-Fisher fixed point 



3 



B, Isotropic O(N) Spin-ice 

When N > 1, the spin takes continuous values, there- 
fore the formalism of the phase transition based on vertex 
operator in the previous section is no longer applicable. 
Also, it is impossible to write down a vertex operator 
with the O(N) spin symmetry. Therefore, we need to 
seek for a different formalism. We consider the following 
Hamiltonian for O(N) spin-ice in addition to the domi- 
nant constraint Eq. |3j 

e = y] JiSi-u,uSi,u + y] j2Sj M sf +UiU . (9) 

i,^j,,a i,a [i^f 

J\ is a Heisenberg coupling between spins along the same 
lattice axis, Ji is a Heisenberg coupling between spins on 
two parallel links across a unit square face. If J\ > 
and Ji < 0, in the ground state spins are antiparal- 
lel along the same axis, but parallel between parallel 
links across a unit square i.e. (S* + „) ~ (— <S° 
is a constant 0(N) vector. Using the CDM terminol- 
ogy, we will call this state the columnar state. If Ji and 
J2 are both positive, in the ground state the spins are 
antiparallel between nearest neighbor links on the same 
axis, as well as between parallel links across a unit square 
{(Sf +/J ) ~ (— lySfy, and we will call this state the stag- 
gered state. Since in the J\ — J 2 model Eq. [9] there is 
no coupling between different axes along different direc- 
tions, for both cases the zero temperature ground state 
of model Eq. [9] has large degeneracy, because the energy 
does not depend on the relative angle between Sf x , S? y 
and Sf z i.e. the ground state manifold has an enlarged 
[O(iV)] 3 symmetry. However, at finite temperature this 
accidental enlarged symmetry of the ground state will be 
broken due to thermal fluctuation. 

In this paper we will focus on the staggered spin or- 
der. Following the magnetic field formalism of the CDM 
mentioned before, in order to describe this system com- 
pactly, we introduce three flavors of 0(N) vector field 
4% = Sf +ll r]i with /i = x,y,z, and now the constraint 
Eq. [3] can be rewritten concisely as 



£v^ = o. 



(10) 



Under the lattice symmetry transformation, 0° trans- 
forms as 

Tfi : fi-tn+1, 4> a v -C, (for all p, u), 

R»,s ■ ft -> -A*, 4>1 -> -4%, C C (for v ^ /i), 

: n <-> v, O cf>l, (for v ± p). (11) 

T M is the translation symmetry along /i direction, R^s 
is the site centered reflection symmetry, and R^ u is the 
reflection along a diagonal direction. 



The staggered spin order corresponds to the uniform 
order of 0®, and all flavors of spin vectors are ordered. 
Therefore presumably (fft are the low energy modes close 
to the transition, and we can write down the following 
symmetry allowed trial field theory for </>^ with softened 
unit length constraint: 



T = £^(-v 2 +r- 7 v 

fj, . a 

+ E-9(E V ^) 2 +^ 



(12) 



When we take g — )• 00, the constraint Eq. \TU\ is effec- 
tively imposed. In Eq. [12] when 7 = 0, the quadratic 
part of the field theory is invariant under O(N) x 0(3) 
transformation, the 0(3) symmetry is a combined flavor- 
orbital rotation symmetry. 7 term will break this sym- 
metry down to the cubic lattice symmetry and O(iV) spin 
symmetry. The flow of 7 comes from the two loop self- 
energy correction diagram (Fig. QJ?), and the RG flow of 
7 will contribute to the RG equation at the order of e 3 , 
which is negligible at the accuracy of our calculation if we 
take e = 4 — d small. Therefore 7 is a constant instead of 
a scaling function in the RG equation, hereafter we will 
always assume 7 is small. 

J- 4 in Eq. [12] includes all the symmetry allowed quartic 
terms of dfl, : 



Fa 



££(^) 2 ] 2 +-£EW) 2 ]£( 

jU a fj,<u a b 



The u and v terms are invariant under an enlarged sym- 
metry [0(./V)] 3 , while the w term breaks this symmetry 
down to one single O(N) symmetry plus lattice symme- 
try. As already mentioned, the ground state manifold 
of model Eq. [9] has the same enlarged [O(iV)] 3 symme- 
try. However, the w term can be induced with thermal 
fluctuation through order-by-disorder mechanism [25| , or 
we can simply turn on such extra bi-quadratic term en- 
ergetically in the J\ — J2 model Eq. [5] By adjusting 
the ratio between u, v and w in Eq. 1131 points with 
various enlarged symmetry can be found. For instance, 
if w + v — 2u, the J- 4 has a O(N) x 0(3) symmetry, 
where the 0(3) is the flavor-orbital combined rotation. 
Just like the Jj — J2 model on the square lattice [25j |. 
the quadratic coupling ^ a <^„0" with /i ^ v as well as 

more complicated quartic terms like [J2 a C^lEt 4>\4>\\ 
break the reflection symmetry of the system, and hence 
are forbidden. 

If we just take the inverse of the Gaussian part of Eq. 
[12] with r > 0, we obtain the following correlation func- 
tion of 0°: 
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Sab 



r + q- 



■(<W 



r + (1 + g)q 



(14) 
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FIG. 1: (a), (b) and (c), momentum-shell loop integrals evalu- 
ated to be j41n(A/A), Bln(A/A) and Cln(A/A) respectively, 
with parameters given by Eq. 1161 (d), the self-energy correc- 
tion that renormalizes 7 in Eq. 1121 this two loop diagram will 
contribute to the renormalization of u, v and w at order e 3 , 
and hence is negligible in our calculation. 



In the limit with g — > 00, the correlation function reads 



g->oo i' ** r + q 

P - x q » q » 

■* flit — u fj,is 2 



(15) 



P^v is a projection matrix that projects a vector to the 
direction perpendicular to its momentum. After Fourier 
transformation, this correlation function gives us the 1 /r 3 
power-law spin correlation of the Coulomb phase. When 
r < 0, the vector 0° is ordered. Interaction T\ will not 
spontaneously generate longitudinal spin wave. For in- 
stance, we can first keep g finite at the beginning, and 
calculate the leading order self-energy correction using 
Fig. Q]l It is straightforward to check that in the result 
after taking the limit g — > 00 the longitudinal wave still 
acquires infinite kinetic energy, and the dressed correla- 
tion function is still fully transverse. 

Now a systematic renormalization group (RG) equa- 
tion can be computed with four parameters u, v, w and 
r, at the critical point r — with the correlation func- 
tion Eq. 1151 In the calculation an e = 4 — d expansion 
is used, and the accuracy is kept to the first order e ex- 
pansion. Based on the spirit of e expansion, all the loop 
integrals should be evaluated at d = 4, and because of 
the flavor-orbital coupling imposed by the constraint Eq. 
[3j we should generalize our system to four dimension, 
and also increase the flavor number to /J, = x, y, z, r. 
Notice that the flavor-mixing correlation in Eq. [15] sig- 
nificantly increases the number of diagrams that we need 
to evaluate. Our RG calculation is based on momentum 
shell integral, since the correlation function Eq. [15] ac- 
quires strong momentum direction dependence, after the 
momentum shell integral the logarithmic correction will 
depend on the flavor of the loop integrals. For instance, 
the loop diagrams in Fig. Q] a, b and c expanded to the 



first order of 7 is evaluated as 
A = 



Ja<| 9 1<A d 4 qD^(q)D^(-q) 



B = 



C 
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(16) 
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7 + 0( 7 2 ). 



For arbitrary N, the full coupled RG equation at the first 
order e expansion reads 



du 
dhd 



dv 
dhd 



dw 
dhd 



dr 
d\nl 



eu - 8A(8 + N)u 2 - 6N(A + 2B)v 2 

6(A + 2B)w 2 - 24(25 + BN)uv 

12(A + 2B)vw - 72Buw, 

ev - 165(4 + N)u 2 - A{2AN + 7BN + AC)v 2 

ACw 2 - 16(2/1 + AN + 2BN + AB)uv 

8(2A + 9B)vw - 16(A + 6B)uw, 

ew - 64Bu 2 - 16Bv 2 

4(2A + IQB + BN + 2C + CN)w 2 

16(5 + 2C)vw - "&2Auw, 

2r - 8(2A + AN + 6B + 'iBN)ur 

12(AN + 3BN)vr-12{A + 3B)wr. (17) 



A similar set of recursion relations of quartic interaction 
terms were computed in a different context in Ref. [261 ] - 
Let us first discuss the solution of this RG equation 
with 7 = 0. Solving this equation at r = with number 
A, B, C given by Eq. 1161 we find eight fixed points, with 
one stable fixed point for N > N c = 70, while for any 
N < N c only instable fixed points are found. The ana- 
lytical expression of the stable fixed point as a function 
of N with N > N c can be straightforwardly obtained by 
solving Eq. [XT] but the result is rather lengthy. Instead, 
we will analyze the solution of Eq. [T7] with an expansion 
of 1/N. For instance the stable fixed point is located at 
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42iV ~ 1323iV 2 
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3e 4870e e s 

W * = 77V- 1323^ + 0( i^ ) ' (18) 

Since now +w* = 2u#, this fixed point has the enlarged 
O(N) x O(cf) symmetry mentioned before. Close to the 
stable fixed point, the three eigenvectors of the RG flow 
have scaling dimensions 



Ai = 
A 2 = 
A 3 = 



-e + 0( 



N 2 ' 
4448e 
567iV " 
24950e 



Q7N 



(19) 



At the stable fixed point Eq. [TH r is the only relevant 
perturbation, and plugging the fixed point values Eq. 1181 
back to the RG equation, we obtain the critical scaling 
dimension 



1 

v 



158e 
W 



(20) 



Since at the ground state all three flavors of spin vectors 
are ordered, in the field theory J-4, v should be smaller 
than 2u, which is well consistent with the stable fixed 
point in Eq. [18] with negative i>* . This fixed point has 
positive Wit , which favors noncollinear alignment between 
spins on different axes. Therefore the transition between 
Coulomb and noncollinear staggered state has a better 
chance to be described by this fixed point. 

Notice that had we included the anisotropic velocity 7 
into account, its leading RG flow will be at order of e 2 , 
and the flow of 7 will contribute to the RG flow of u, v 
and w at order of e 3 , therefore it is justified to take 7 
a constant in our calculation as long as we keep e small 
enough. When 7 is nonzero but small, we can solve the 
RG equation with A, B and C given by Eq. [T51 and 
the RG flows will only change quantitatively, although 
the 0(N) x 0(3) symmetry of the stable fixed point is 
broken by 7. Expanded to the first order of 7, the scaling 
dimensions of the three eigenvectors of the RG equation 
at the stable fixed point become Ai = — e, A2 = — e + 



4448e , 2849936^7 



A, 



249506 30088528e 7 nTlr | t >.„ 
27921915W ' a " 



567N 1 3988845^ ' ^ J — c ^ 567N 

scaling dimension of r becomes [r] = ~ = 2 — e + — 

103267 v 
X715N ■ 

If we took the limit TV — » 00 in the physical system, 
we only need to keep the terms linear with N in Eq. 1171 
and now the equation becomes precise even with e = 1. 
In this case the RG flow of w is decoupled from u and v, 
hence four of the eight fixed points have = 0, and all 
the others have u>* = 3/(7N). The RG flow diagram for 
u and v with w = w* = 3/(7N) in the large- N limit is 
depicted in Fig. [5Ji. 

As we promised in the beginning of this paper, we 
should discuss the applicability of the constrained GL for- 
malism discussed in this paper to the CDM-3, which cor- 
responds to the case with N = 1. In our GL formalism, 



in the ordered phase, the power law spin-spin correlation 
still persists if the long range correlation is subtracted. 
For instance, the fluctuation <50 M = (f)^ — ((f>^) is still sub- 
ject to the constraint ^2 = 0, therefore although 
the fluctuation is gapped, it still leads to the 1/r 3 power- 
law correlation. But in CDM, the ordered phase only 
has short range connected dimer correlation on top of 
the long range order, which can be checked with a low 
temperature expansion of CDM ■ Like what was dis- 
cussed in the introduction, the key property of the case 
for N — 1 is that, the spins only take discrete value ±1, 
the vertex operator like Eq. [5] in the dual field theory in 
terms of vector potential can drive the system to a phase 
with short range connected correlation through a Higgs 
transition. The effect of the vertex operator was missing 
in our GL formalism. 

We can solve Eq. [T7] with N = 1, where v and w 
terms are identical. In this case in addition to the trivial 
Gaussian fixed point, there is only one other fixed point 
at = 2u* = e/34 with 0(3) flavor-space combined 
rotation symmetry, which is the same fixed point as the 
ferromagnetic transition with dipolar interaction 2|| 2t| . 
In 3d space, the dipolar interaction also projects a spin 
wave to its transverse direction. The dipolar fixed point 
is instable against the 0(3) to cubic symmetry breaking, 
therefore when N = 1 our first order e expansion predicts 
a first order transition. Based on the discussion in the 
previous paragraph, one possible scenario for the CDM- 
3 with staggered ground state is that, if we lower the 
temperature from the Coulomb phase, after the first or- 
der transition of 4> a , there has to be another "Higgs" like 
phase transition that destroys the power-law connected 
correlation. Or there can be one single strong first order 
transition that connects the Coulomb phase and stag- 
gered dimer crystal directly. 



C, O(N) Spin-Ice with Easy Plane Anisotropy 

Model Eq. [9] is invariant under cubic lattice symmetry, 
and we can certainly turn on various anisotropies to this 
model, like what was studied in the CDM-1 model 
For instance, let us modify the model Eq. [9] slightly: 



E 



/j ?j J"2,nuSf ^S°; +u l l2l) 



If Jl,x = Jl.y > Jl.z > 0, and J2,xy = J2,xz = J2,yx = 

J"2,yz > J2.ZX = J2,zy > 0, the O(N) spin vectors in 
the xy plane, Sf x and Sf y have a stronger tendency to 
order compared with Sf z . Therefore when we lower the 
temperature from the high temperature algebraic phase, 
the O(N) vectors in xy plane are expected to order first 
at critical temperature T c \. In the field theory close to 
T c i, the anisotropy can be described by an extra mass 
term for ^2 a rn 2 ((j)^)' 2 in the free energy, which is clearly 
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FIG. 2: (a), the lattice structure of the O(N) spin- ice model. 
(6), one of the plaquette configurations of CPM-2, with every 
site shared by two plaquettes. (c), lattice structure for the 
O(N) spin-plaquette model, each unit square face of the cubic 
lattice is occupied by an O(N) vector, and the constraint is 
that the sum of all the O(N) vectors on the twelve square 
faces shared by site i is zero, (d), the RG flow of Eq. 1171 in 
the large-iV limit, C is the stable fixed point, with enlarged 
O(N) x 0(d) symmetry. 



a relevant perturbation at the critical point r = 0: 

+ E*C 



(E 

A* 



V 



(22) 



In the equation above we have taken 7 = for simplic- 
ity. To calculate the RG equations for J4, we still need 
to increase the dimension and flavor number to four with 
fj, = x, y, z, r, and the anisotropy of the generalized 
system will prefer the O(N) vectors to order on three of 
the four axes. At the critical point r = 0, due to the rele- 
vance of the extra mass term, we can safely take m — > 00, 
and the correlation function between </>° reads 



lim 

Q His 



;>M(-9)>~^ 



Suv 



..2 



= 0, (i or v = t. 



Using the correlation function Eq. 
at the critical point reads 



C + ^ + 3t! 



(23) 

the RG equation 



du 
dhTl 



eu - 8A(8 + N)u 2 - AN{A + B)v 2 



dv 
dhd 



dw 
dhd 



dr 
dhd 



4(A + B)w 2 - 45(8 + 4N)uv 
8(A + B)vw - A8Buw, 

ev - 165(4 + N)u 2 - A(AN + 3BN + AC)v 2 
ACw 2 - 16(2,4 + AN + BN + 2B)uv 
8(A + 5B)vw - 16(A + 3B)uw, 



ew - 6ABu 2 - 16Bv 2 



- 4(A + 6B + BN + 2C + CN)w 2 

- 16(5 + 2C)vw - 32Auw, 

= 2r -8{2A + AN + AB + 2BN)ur 

- 8{AN + 2BN)vr — 8(A + 2B)wr. (24) 



Now A = 8/15, B = 1/15, C = 2/5, which is different 
from the isotropic case, due to the different form of the 
correlation functions. Solving this equation at r = 0, 
we again find a stable fixed point for large enough N. 
Expanded to the order of e/N 2 , the stable fixed point is 
located at 

27e 1431e 



112AT 49CW 2 
3e 879e 
~56N ~ 490N 2 
15e 1983e e 
28N ~ A90N 2 + ^TV^' 



(25) 



Again, since + w* = 2m*, this fixed point has an en- 
larged O(N) x 0(d — 1) symmetry. Close to the stable 
fixed point, the three eigenvectors of the RG flow have 
scaling dimensions 



942e n . e , 
= - £+ 25A7 + 0( ^ 
1104e e . 



(26) 



A2 is the largest scaling dimension, and according to Eq. 
[24]the critical N is N c ~ 59, below which the stable fixed 
point disappears. Plug the fixed point value Eq. 1251 back 
to the last RG equation in Eq. fJH we obtain the scaling 
dimension of r at the fixed point: 



1 

v 



138e 
~7N~ 



N 2 ' 



(27) 



If we keep lowering the temperature after the order of 
<j>% and <j>y, then eventually 4>t will also order at temper- 
ature T c2 < Tax. After the order of #J and the O(N) 
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symmetry is broken down to its subgroup. Let us first 
assume w < Q in .F4 i.e. the three flavors of 0(N) vec- 
tors are collinear with each other in the low temperature 
ordered phase. Let us assume that in the intermediate 
phase the expectation value (0") ~ (cj)y) ~ (1,0, •••0), 
now the O(N) symmetry is broken down to 0(N— 1) sym- 
metry generated by Lie algebra L a b with a, b = 2, • • • N. 
This symmetry breaking implies that the degeneracy of 
the N components of <f> z will be lifted at T c2 . Because 
w < 0, the nonzero expectation value of <\>\ and <fi y 
will prefer <p\ to order next at lower temperature, which 
is essentially an Ising transition with order parameter 
(f>g. From now on we will denote 
couple to the gapped fluctuations < 



as 



(j) z will 
(<f>l) and 

<j)y = <j>y — (<j) y ) through the constraint, and the entire 
low energy field theory at T c2 reads 



E ^(-V 2 + m 2 )</. p 



.9(E V ^^) 2 + 0(^ 



(28) 



After taking the limit g — ¥ oo, the critical correlation 
function of <A, reads 



lim D zz (q) 

g->oo 



1 



(29) 



Therefore this transition at T c2 is effectively a z = 2 tran- 
sition, with scaling dimension [q z ] — 2\q x ] — 2[q v ] = 2. 
Now the total effective dimension is 4, and the (0 Z ) 4 
term is a marginally irrelevant perturbation, therefore 
this transition is a mean field transition with logarithmic 
corrections. Notice that the Goldstone modes after the 
O(N) to 0(N — 1) symmetry breaking are harmless to 
this transition. The Goldstone mode can be described 
by 8(j>% ~ 8<j)y ~ (0, 7r 2 , • • • tt n ) that forms a vector rep- 
resentation of 0(N — 1), and in order to guarantee the 
gaplessness of the Goldstone modes, close to T c2 the fol- 
lowing coupling between <j) z and 7? is the lowest order 
coupling that is allowed: 



J 7 no 



Idsto 



E« 



^(V M 7?) 2 



(30) 



This term only generates irrelevant perturbations for <f) z 
at the transition. Couplings like 0^(7t) 2 is forbidden due 
to its ability to renormalize the mass of Goldstone mode 

7?. 

If w > in then the three flavors of O(N) vec- 
tors are perpendicular to each other in the low tempera- 
ture ordered phase. Let us assume that in the interme- 
diate phase between T c \ and T c2 the expectation values 
(4>%) ~ (1,0, •••0) and - (0, 1,0 •••0). Now the 

O(N) symmetry is broken down to the 0(N — 2) sym- 
metry generated by F a ;, with a, b = 3, • • ■ N. Due to the 
presence of w term in J-4, at T c2 the order parameters 



should be § a z with a = 3 ■ • • N, which forms a vector rep- 
resentation of 0(N — 2). In the intermediate phase there 
are in total 2N — 3 Goldstone modes, they are 7r" = 
(0, 0, 7r|, • • • ir^f ) that corresponds to the O(JV) Lie alge- 
bra elements Ti a with a = 3, • • • N; 7r™ = (0, 0, 7r^ • • • ) 
that correspond to the the O(N) Lie algebra elements T 2a 
with a = 3, • • • N; plus a rotation mode between <j>\. and 
</> 2 , denoted as F12. Again the order parameter </3" will 
couple to these Goldstone modes through the constraint. 
The field theory close to T c2 reads 

N N 
fj.=x,y a=3 



a— 3 



N 



+ E -9( V ^ + V *< + V ^0 2 + (31) 

a— 3 

Notice that the interaction between Goldstone modes 7r° 
have to be irrelevant to guarantee their gaplessness. The 
universality class of this transition at r = can be cal- 
culated using the correlation function of <fc® after taking 
the limit g — > 00: 



lim DfM 



ill 



8 a b 



(32) 



The first order e expansion leads to the following scaling 
dimension and critical exponent: 



iV — 1 
N + 5 



(33) 



which is identical to the 3d 0(JV-2) Wilson-Fisher (WF) 
fixed point, although higher order expansions may de- 
viate from the WF fixed point. Therefore our formal- 
ism implies that with large enough N, the easy plane 
anisotropy will split the transition discussed in the pre- 
vious subsection to two second order transitions. 



D, O(N) Spin-Ice with Easy Axis Anisotropy 

In model Eq. [2TJ if we make the following choice of 
parameters: J^ z > Ji jX = J\ tV > 0, J 2 ^ zx = J 2 , zy > 
■h,xy = -h,xz = J2, V x = J%,yz, then when we lower the 
temperature from the algebraic phase, the O(N) vectors 
along the z axis will order first at temperature T c \ . Sim- 
ilar to the previous section, this easy axis anisotropy can 
be described by an extra mass gap m for both <f)% and (j> y 
modes in the field theory Eq. [T2J and (f> z becomes the 
only order parameter at low energy. However, now we 
can no longer take the limit m — > 00, because with this 
limit all the correlation functions will vanish. If we keep 
m finite, the correlation function of <$>\ takes a similar 
form as Eq. \M 



lim DfM 



Sab 



nv 



it 



(34) 
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therefore this transition is again an effective a z = 2 
mean field transition. A similar situation with easy axis 
anisotropy has been studied in Ref. [3(|, although there 
the easy axis was along the diagonal direction. 

If the temperature is lowered even more from T c \ , then 
at T C 2 the spin vectors within xy planes will order next. 
Again the nature of this transition would depend on the 
sign of w in Fi. If w < i.e. the three flavors of O(N) 
vectors order collinearly, then if 0" orders along direction 
(0°) ~ (1,0- --0), at temperature T c2 4> x and 4> y will 
become the critical modes. 4> x and <pz are coupled to the 
gapped fluctuation <p z = <j>\ — (4>\), and the low energy 
field theory describing cj> x and <p y is identical to Eq. [22] 
with N = 1. By solving the RG equation Eq. [24] with 
N = 1, no stable fixed point is found. 

If w > 0, the critical modes at T c2 is <p® and <j>y 
with a — 2, • • • N, which form vector representations of 
0(iV— 1). These two critical modes are coupled to gapless 



Goldstone mode 7r" = (0, tt z 



with field theory 



N 



N 



fi—x,y a— 2 a— 2 

N 

+ ^.9(V z vr 2 a + V x 4>% + V,^) 2 + J- 4 (C, 02X35) 



a— 3 



Now T± only involves two flavors of 0(iV — I) vectors 
and 4>y. Number K is in general not 1 because the 
now there is no symmetry that can transform <p x , <fiy and 
7r°. Again, the renormalization of K will be an order e 2 
effect, therefore we can take K as a constant. The RG 
equation for J-4 can be calculated in the same manner as 
we did before. If we take K = 1, then the RG equation 
for .F4 takes the same form as Eq. [M]but with A = 5/8, 
B = 1/24, C = 13/24. With large enough N there is 
still a stable fixed point, and the answer is qualitatively 
unchanged when K deviates from 1 slightly. 



Ill, 0(N) SPIN-PLAQUETTE MODEL 
A, Phase Transition with General N 

Now we switch the gear to the less well studied spin- 
plaquctte model, where the physical quantities are de- 
fined on the unit square faces of the cubic lattice instead 
of the links. In Ref. 1|| TtJ, a model of this type was 
studied, with every site of the cubic lattice connects to 
precisely one filled plaquette. Physically the filled pla- 
quette is a SU(4) singlet formed with four SU(4) funda- 
mental fermions on four corners of the plaquette. This 
model is denoted as CPM-1. In this section we will study 
an O(N) generalization of this model. We define an O(N) 
spin vector Sf ^ v on each unit face of the cubic lattice, 



and impose the constraint 

^ ] ^i,[ii> ^i— p,,p,v ^i—v.p.v ^i— p,— v^v 0.(36) 



In addition to this constraint, we consider the following 
Hamiltonian: 

E = E ^^i,itv^i+n,nt> + J2Sf^ v S°; +p ^ v ,p ^ p, v. (37) 



When Ji > 0, J2 > 0, the ground state of this 
Hamiltonian has staggered order with nonzero ((f>f „„) ~ 
(— iySf u with p, v = +x,+y,+z. is a symmetric 
tensor with p ^ u, which has only three independent fla- 
vors. Under discrete cubic lattice symmetry, tensor field 
^„ transforms as 



T„ : i^i + p, § a vp -)• -<f>" p , (for all fx, v, p), 



R,. 



p -> -p., 4>l v -»• -tf^,, ft -> 4> a vp , (for u,p + p), 



(for v± p.). 



(38) 



Just like the spin-ice model Eq. [H] the ground state of 
model Eq. [37] also has enlarged symmetry [O(iV)] 3 i.e. 
the spins on different planes will order independently. 
Since <$> a is subject to the constraint Eq. [5] we can start 
with the following low energy field theory close to the 
transition: 

f = E ^(^ v2 +^+Eff(E v ^ v ^) 2 

+ J" 4 , (39) 
and J-4 is all the quartic terms allowed by symmetry: 

* = -EE(^) 2 ] 2 +-EE(^) 2 hE«p) 2 ] 

b 

(40) 



p<u a 

+ ^EE<M4E^ 

v^p a b 



u<p a 

6 1 



Under the limit g — > oo, the constraint Eq. [6]is effectively 
imposed, and the correlation function with r > is: 



r + q 



2 1 ^fwpa • 



QxQy + lllz + 9 2 9 2 ' 



(41) 



Notice that in Eq. [39] some extra flavor-orbital coupling 
quadratic terms are allowed by symmetry, such as 

^2~E V ^ a y^v^ + Vy€ x VA a xy + V Z C,V^142) 



Just like the 7 term in Eq. [T^l this term will not gain 
any renormalization at the one-loop calculation i.e. its 
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RG flow will only affect the RG equation of u, v and w 
at the order of e 3 . A small perturbation of this term will 
not change the RG flow qualitatively, we will take this 
term to be zero hereafter for simplicity. 

To calculate the RG equation at the critical point 
r = 0, in principle we need to generalize the model to 
4d. However, there is no ideal way to make this general- 
ization. For instance, if we define O(N) vectors on the 2d 
faces of a 4d lattice, and ju takes x, y, z, r, at 4d there 
are two extra quartic terms in addition to T± in Eq. 1401 



4(7. 



a, 



E ^i/^E H 

a b 

E^v«rlE^™] + --- (43) 

a b 



These two extra terms will make the RG equation much 
more complicated than the actual 3d case. For instance, 
since these two terms are allowed by symmetry, they will 
be generated under RG flows even if we take them to be 
zero at the beginning. Therefore in this section we will 
just evaluate the loop integrals at 3d. 

The RG equation takes exactly the same form as Eq. 
[Ml while now A = 5/9, B = 1/18, C = 7/18. With large 
enough N, there is a stable fixed point located at 



lie 

487V 
e 



29579e 
10368iV 2 
19039e 
24iV ~ 10368iV 2 
9e 22171e 



!(>* 



167V 5184iV 2 



0( 



N3> 



(44) 



Close to the stable fixed point, the three eigenvectors of 
the RG flow have scaling dimensions 



Ai = -e 



A 2 = -£ 



6.80e 



N 
37.5e 

N 



0(-' 



A, = - e+ ^ +0 (_). (45) 

Plug the fixed point value Eq. [55] back to the last RG 
equation in Eq. 1241 we obtain the scaling dimension of r 
at the fixed point: 



W = I = 2 



2113e 



108N ' ~ y N 2/ ^ 

UN — 1, no stable fixed point is found. The case N = 
1 is equivalent to CPM-6, which is much less studied 
compared with the CDM-m. We will discuss this model 
in the next subsection. 



B, Duality and N = 1 

Besides taking the limit g — !> oo in the Gaussian field 
theory Eq. [321 the correlation function in Eq. 2JJ can be 



obtained by other means. Let us solve the constraint Eq. 
[6] with softened unit length constraint of <p" by defining 
the three flavors of O(N) height field on the cubic 
center i of the lattice: 



K) 



<t>1x = VAK-K)- 



(47) 



Then height field plays the same role as vector poten- 
tial in the usual dimer model. Now in terms of the 
height field, the algebraic phase can be described by the 
following Gaussian field theory 



T = 



E- 

^ 2 



[v z {K-h a y )f + [v y {K-K) 



(48) 



This Gaussian field theory takes a similar form as the 
Khaliullin model describing the orbital degrees of free- 
dom [31(, after taking the spin- wave expansion. Diago- 
nalizing this Gaussian field theory, we obtain two eigen- 
modes describing the fluctuations of the height vector 
field: 



2 2 

Qxly 



(49) 



Interestingly uji vanishes along each coordinate axis in 
the momentum space. This height field theory Eq. 1481 is 
invariant under the following symmetry transformation: 

K + h{x) + ip{x,y,z), 



h a 



h a 

fly 



h(y) + <p{x,y,z) 



K -> K + f 3 (z) + cp(x,y,z). 



(50) 



f a are functions of only one of the coordinates, this type 
of quasi-local symmetry comes from the definition of the 
height field Eq. \$7\ and hence does not depend on the 
detailed form of the Hamiltonian. This quasi-local sym- 
metry is responsible for the line of nodes of the eigen- 
modes in Eq. 03 Using the Gaussian free energy, and 
the height representation of </>^„ in Eq. 1471 we can repro- 
duce the correlation function Eq. ST] 

These dualities are particularly useful for N = 1, which 



is equivalent to CPM-6. Since in this case 



takes 



discrete values, then are also discrete. Therefore 
in the dual field theory Eq. [48] for N = 1 we should 
also consider the vertex operators like we introduced for 
the CDM, in section II-A. For convenience, let us make 
the following standard modification of our description 
of the system: we will allow cj)^ to take all the half- 
integer values, and turn on classical Hamiltonian on the 
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lattice E — Y^i a^u U(<f>i,fj,i>) 2 , on top of the constraint 
Eu#w ^ttVv<j>nv = 0. When U is large, effectively on 
every unit face 4>i t uv can only take two values, which is 
the same as the CPM-6. In this way, even though we 
increased the total configurations of cj)^, the low energy 
configurations of <^„ is still identical to the CPM-6. Now 
the dual field theory of the algebraic phase reads 



[Vy{h z -h x )f 



(51) 



In order to make sure (j)^ take half-integer values, based 
on the definite Eq. [47] — h v will take site-dependent 
integer or half-integer values on the 3d dual cubic lattice: 



h x — h y G Z ■ 



l-(-l) 



h y — h z G Z ■ 



h z -h x G Z 



I- {-If 



l-(-l) 



i v +iz 



(52) 



Then in the dual theory the lowest order vertex operator 
that do not have spatial oscillation is 



T v = 2J -acos[47r(/i Al - h v )}. 



(53) 



It was shown in Ref. [17| that this vertex operator has 
directional dependent algebraic correlation in the alge- 
braic phase, with scaling dimension proportional to K in 
Eq. [5TJ For instance, let us denote cxys[Air{h x — h y )] as 
J~v, X y, then due to the quasi-local symmetry in Eq. 150] 
two T v ^ X y operators can only have nonzero correlation 
when they are on the same z axis: 



(^(0,0,0)^(0,0,2)) 
4(2tt) 2 f d 3 k ikl 



exp{- 
1 

z 16tt/K ' 



kl)e 



ik z z 



} 



(54) 



If the scaling dimension 8tt/K is greater than 1 i.e. 
K > K c i = 8ir, the vertex operator Eq. [53] becomes 
relevant in the algebraic phase. The relevance of this 
vertex operator can be manifested by directly calculat- 
ing the partition function with expansion of the vertex 
operator. The second order perturbation involves the in- 
tegral of the correlation function Eq. [54] and if K > K c \ , 
the integral diverges in the infrared limit. 

the proliferation of the vertex operator will drive a 
Kosterlitz-Thouless (KT) transition, across which the al- 
gebraic correlation disappears, and the correlation length 

of the vertex operators diverges as £ ~ expL I K _^ K ]• 



This KT transition in 3d space is due to the special 
quasilocal symmetry in Eq. [50] which grants the 3d sys- 
tem a 2d like symmetry for each flavor of h^. The KT 
transition and dimensional reduction behavior in 3d or 
2+ Id were also discussed in another type of U(l) ro- 



tor systems with similar quasilocal symmetries [32h34 | . 
What is different here is that, this dimensional reduc- 
tion behavior in our system is inherited from the generic 
quasi-local symmetry in the definition of height field Eq. 
l47l therefore dimensional reduction is robust. 

There is another useful duality of Eq. 25] Consider 
the following Gaussian field theory: 



N 1 

a=l 



a\2 



(V y V z a y + (V z V^ a )1(55) 



Now let us introduce the new field (ff^v through Hubbard- 
Stratonovich transformation: 



K 



(56) 



After integrating out the field 9 a , the partition function 
of the system becomes 



Z = J D^jCZ^^^e-^ 3 ^^ 



"(57) 



The delta function in this partition imposes the same 
constraint Eq. [5] The field theory Eq. [SSlis also invariant 
under quasi-local symmetry transformation: 



na . rta . a 



+ 9l{x) + g a 2 (y) + g a 3 (z). 



(58) 



g° are arbitrary functions of one of the three coordinates. 
The field theories Eq. [S5] and Eq. 05] are dual to each 
other, with the same parameter K. 



IV, DEFECTS OF CONSTRAINTS 
A, Affects on the Algebraic Phase 

So far all the constraints have been perfectly imposed, 
and in this section we will consider the case with slightly 
softened constraint |35[ i.e. we allow the existence of the 
point defect that violates the constraint. In the algebraic 
phase, for all the O(TV) models considered in this work, if 
g is large but finite in Eq. [T^] HH EE the power-law cor- 
relation immediately crossovers back to short range cor- 
relation for large enough distance. For instance, take the 
correlation function Eq. [M] we can see that the length 
scale for this crossover is / ~ \J gjr for O(N) spin-ice. 

The Gaussian theory evaluation of the defects is based 
on the assumption that the defect can take continuous 



values, therefore the "charge" of the defect J^fj, can 
be infinitesimal. Therefore the Gaussian theory is no 
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longer applicable for N = 1, where the defect charge is 
always discrete. To evaluate the defect in this case, one 
needs to go to the other side of the duality. For instance, 
the CDM-m is dual to the Villain form of U(l) rotor 
model with partition function 



Z = 



exp[ 



Notice that there is an imaginary term in the partition 
function due to the nonzero average filling of the dimer 
density, which is similar to the Berry phase of quantum 
Bose rotor model with fractional boson filling [36|. To 
show this duality explicitly, we still introduce the field 
<pi ttl through Hubbard-Stratonovich transformation 



Z = J DOD^^MzZ 

■ Tfl 

- 27ril il/l (& l/1 + ("^V )]• 



K 



2 flu ~ * V /A/A 



(60) 



After integrating out the 0j, and summing over li tfjt , the 
dual theory reads 



K 



) 2 
1 1 



(61) 



1 ■/' 



and (— iy<f)i jfl € Z + m/6. Then any low energy configu- 
ration of cf>i^ is equivalent to a CDM-m configuration. 

If we are in the superfluid phase of the rotor model, 
which is dual to the algebraic phase of the CDM-m, we 
can expand the free energy of the rotor model Eq. [55] at 
0, and the field theory of the superfluid phase is 
) 2 . The defects will be taken into 



simply E = 2f(V m 0) 2 
account by the vertex operator in the field theory: 

^ = E^( v ^) 2 -« cos W> 



(62) 



a ~ exp(— C\g) is the fugacity of the defect, and the 
partition function with expansion of a is equivalent to 
a classical Coulomb gas, which is equivalent to the par- 
tition function of defects. In 3d space since 9 has long 
range order, a is a very relevant perturbation i.e. the 
presence of defects will immediately destroy the algebraic 
phase with infinitesimal fugacity. Another way to show 
the relevance of the defect is to compare the energy and 
entropy of an isolated defect. If the system size is L 3 , 
then an isolated defect in the algebraic phase costs en- 
ergy E ~ J d 3 kj% ~ gA+j; which is finite in the infrared 
limit. While the entropy of the defect scales as S ~ InL, 
therefore the entropy always dominates energy i.e. the 
defects always proliferate. 

The situation is very different for the CPM-m. Just 
like the duality of the CDM-m, the CPM-m is dual to 



the following rotor model: 



Z = 



f £>0£exp[ - 



E ^(V AI ^-27T^) 2 



. m 



iTri-kM-iy 



(63) 



Again, in the algebraic phase of the CPM-m, we can 



expand the free energy at 



0, and evaluate the 



(59) relevance of defects in the rotor field theory: 



— a cos(0). 



(64) 



Now because of the quasi-local symmetry Eq. [551 the 
vertex operator J- v = —a cos(#) has no nonzero correla- 
tion spatially, which seemingly implies that the defect is 
irrelevant. This effect can again be shown by evaluating 
the energy and entropy of the defect. An isolated defect 
costs energy 



E 



d 3 k 



1 



£-2 i>2 



L + clni, 



(65) 



which always dominates the entropy that scales as In L 
i.e. the defects are always suppressed by energy. 

However, although nonzero correlation between !F V is 
forbidden by symmetry, the correlation between defect- 
dipole operator Ed,^ ~ — a'cos(V M #) = — a' cos(#i — 
0i+^) can be nonzero. For instance, Ed, z can have 
nonzero correlation within the entire xy plane: 



(E d , z (0,0,0)E d}Z {x,y,0)) 



exp{— K 
1 



d 3 k 



fc2 ^ik^x+ikyy 



h2h2 i 1.21,2 i h 2 h 2 



} 



{x 2 + y 2 ) K /(^) 



(66) 



When K < K C 2 — 8tt, the vertex dipole operator E d ^ 
is relevant, and the algebraic phase disappears. Sup- 
pose a' ~ a 2 ~ exp(— C^g) is the fugacity of the de- 
fect dipole, the crossover length scale beyond which the 
dipole becomes important is £ ~ (a') 47T ^ Kc2 ~ K \ Again 
we can compare the energy of a defect dipole and its 
entropy. The energy of a defect dipole scales as E ~ 
/ d 3 kk 2 l(k 2 k 2 + k 2 k 2 + k 2 k 2 ) ~ InL, which is compara- 
ble with the entropy. So by tuning K in Eq. I64( there 
will be a KT transition as a result of the competition 
between entropy and energy. 



B, Effects on the Transition 

Close to the transition r — 0, we can use the correlation 
function Eq. [14] to compute the RG equation, g may 
flow under RG eventually, but it remains a constant at 
the first order e expansion, and hence we will just take 
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it a constant. For the isotropic O(N) spin-ice, the RG 
equation takes the same form as Eq. [TT1 but now we need 
to reevaluate the loop integrals: 



A = 
B = 
C = 



1 



1 



(I + 3) 2 4(l+ 5 ) ! 



24(l + 5 )2' 
13 1 



24 24(l + 3) 2 12(1 + 5)' 



(67) 



Solve the equation Eq. [17] with the new parameters 
Eq. [67] we can see that when g is large, the solutions 
are qualitatively unchanged from Eq. [25] While when 
g < g c ~ 4.1, stable fixed points are found with N = 1, 
which corresponds to CDM-3 with softened constraint. 
For instance, with N = 1, g = 1, the stable fixed point 
is located at = 53 ~^ , + = 7 +^p . Now in 
the low temperature phase the connected correlation on 
top of the long range order is also short-ranged after the 
transition described above. 

With finite g, the O(N) spin-ice with easy-plane axis 
can be studied in a similar way as last paragraph. As 
was mentioned in section II. A, we need to generalize the 
system to 4d with fi = x, y, z, r, and anisotropy of the 
Hamiltonian prefers the O(N) spin vectors on x, y, z 
axes to order first at T c \. At the critical point T c \ with 
r = 0, the 4d correlation function between <$> a reads 



m— >oo,r— ^ 



gq^q v 



q 2 ^" -Z 2 + (l+3)q 2 



q = 



fl2 



ql 



lim D 



ah 
in/ 



0, \i or v = 0. 



(68) 



The RG equation in this case takes the same form as Eq. 
|2"41 with different A, B and C . Expanded to the first 
order of 1/g, these parameters read: 



.4 
B 
C 



15 153 

1 4 
15 ~~ 153 

2 16 

5 + 15^" 



g 2 

-o{- 2 ). 
g 



(69) 



When 3 < 3 C ~ 4.4, stable fixed points are found with 
N = 1, which corresponds to CDM-3 with softened 
constraint and easy plane anisotropy. With easy-axis 
anisotropy and finite g, the effective z = 2 physics in 
Eq. [34] is absent, and the universality class is expected 
to crossover back to a 3d Wilson-Fisher transition. 

In the O(N) spin-plaquette model, since the g term 
has higher derivatives compared with the Gaussian part 
of the field theory Eq. 1391 then g becomes irrelevant when 
its initial value is finite. Since at the low energy effective 



field theory, there is no flavor-orbital mixing interaction, 
we can take a more concise notation: 0* = (f)y Z , <f)y — 0™ x , 
4>1 = (f>® y . At the critical point r — the O(N) spin- 
plaquette model is described by the field theory: 



■ft, 



(70) 



with T\ given by Eq. [T3] Now the coupled RG equation 
for it, v and w reads 



du 
dhd 

dv 
dhd 



dw 
dhd 

dr 
dhd 



eu - 8(8 + N)u 2 - 4Nv 2 - Aw 2 - 8tw, 

ev - 4(4 + N)v 2 - Aw 2 

16(2 + N)uv - 2Avw - l&uw, 

ew - 4(3 + N)w 2 - I6vw - 32uw, 

2r - 8(2 + N)ur - 8Nvr - 8wr. (71) 



For large enough N, there are in total eight fixed points. 
For instance, the fixed point w* — 0, = 2u* 
< (12N) + - ■ ■ has enlarged 0(3^) symmetry, which mixes 
the flavor and spin symmetry. Expanded to the order of 
e/N 2 , the only stable fixed point with r = is located at 

_ e 5e e 
U * ~ 87V~4iV2~ + 
3e 



w* = — 



4iV 2 1 °( N 2)> 



e 
47V 



7e 



0(- 



(72) 



Close to the stable fixed points, the three eigenvectors of 
the RG flow have scaling dimensions 

e , 



Ai 



-e + 0( 
4e 



N 2 ' 



20e _ e . 



The scaling dimension of r at this fixed point is 



W = i = 2-e. 



12e 



°(^' £2 )- 



(73) 



(74) 



If we take N — 1 , the field theory Eq. [70] becomes 
the standard 0(3) transition with cubic anisotropy, and 
various numerical methods have confirmed that the 3d 
0(3) symmetric fixed point is stable [13], with critical 
exponent [ r ] = l/v = 2- 5e/ll + 0(e 2 ) (see Ref. (H). 



V, SUMMARIES AND DISCUSSIONS 

In this work we studied the classical phase transition 
between the algebraic phase and low temperature spin or- 



13 



dered phase in several different types of 0(A) spin mod- 
els with local geometric constraint. Effects of softened 
constraints are also considered in all the models. Sys- 
tematic RG calculations are applied to all of the cases, 
and solutions at the first order e expansion were obtained 
as precisely as we could. However, higher order e expan- 
sions, as well as direct Monte Carlo simulations of the 
lattice models are indeed demanded in order to confirm 
our results at a more quantitative level. 

So far we have been focusing on the staggered spin 
order in all cases we studied, which has a straightfor- 
ward order parameter description. However, there are 
another large class of spin orders of constrained systems 
that seem to involve more complicated order parameter 
descriptions. For instance, in the columnar order of spin- 
ice, <^ does not have nonzero expectation values, there- 
fore we need to develop another formalism for this case. 
The columnar order is equivalent to order of at mo- 
mentum Qi = (0,7r, 7r), <f>y at momentum Qi — (tt, 0, tt) 
and <f>" at momentum Q\ = (tt, tt, 0). Therefore pre- 
sumably we could describe this transition with conden- 
sation of 4>1j, at all three wave- vectors. However, just like 
the Coulomb-columnar transition in CDM-1 discussed in 
section II-A, some topological configuration of these or- 
der parameters may be forbidden, which potentially can 
change the universality class completely. 

Another interesting subject is to generalize our formal- 
ism to the quantum case. For instance, it is well-known 
that the quantum 0(A) rotor model and spin models can 
be described by nonlinear sigma model in the infrared 
limit. Suppose we impose a local constraint on the order 
parameters of the quantum rotor model or spin model, 
it is possible that we can use field theories similar to Eq. 
[T2l [22l to describe the quantum phase transitions in the 
constrained Nonlinear sigma model: 

a fi—x,y--- 

+ £fl £ v ^) 2 + --- (75) 

with the limit j3 — > oo. Notice that flavor /i takes only 
spatial coordinates. 
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